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Abstract 

The Tokyo group has performed the first underground dark matter search experi- 
ment from 2001 through 2002 at Kamioka Observatory (2700m. w.e). The detector is 
eight LiF bolometers with a total mass of 168g and aims for the direct detection of 
WIMPs via spin-dependent interaction. With an exposure of 4.1 kg-days, we derived 
the limits in the Op — On (WIMP-nucleon couplings) plane and excluded a large part 
of the parameter space allowed by the UKDMC experiment. 
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1 Introduction 



Weakly Interacting Massive Particles (WIMPs) are thought to be the most 
promising dark matter candidates and many experiments have been per- 
formed for direct detection of WIMPs. The DAMA experiment reported hav- 
ing detected WIMPs by an annual modulation signal with their 10kg Nal 
detectors[l]. The reported WIMP parameter region is Mwimp = {52^l^)GeV 
and (T^iMP-p = (7.2lo:9)-10"^pb. On the other hand, the CDMS CoUaboration[2] 
and EDELWEISS Collaboration [3] reported having excluded most of the pa- 
rameter space claimed by the DAMA experiment within the standard su- 
persymmetry(scalar couplings) and the halo models. These experiments aim 
to detect WIMPs via the spin-independent interactions. In 2001, the DAMA 
group published a paper in which they interpreted the annual modulation sig- 
nal in a spin-independent and spin-dependent mixed framework[4]. According 
to Ref. [4] 

) ^wiMP-p could be as high as 0.6pb for 30 ~ 80GeV if we assume cer- 
tain sets of parameters. These results suggest that WIMPs could be detected 
via spin-dependent interaction and the dark matter search experiments should 
be performed both via spin-independent and spin- dependent interactions. 

We, the Tokyo group, had developed a lithium fluoride (LiF) bolometer to de- 
tect WIMPs via spin-dependent interaction in the surface laboratory[5] then 
performed pilot runs at the shallow depth site [6]. We have improved our LiF 
bolometer since we installed the detector system in 2000 at Kamioka Obser- 
vatory (2700m.w.e) and recently we made a dark matter search measurement 
from 2001 through 2002. 

In this paper, we report our first results from the dark matter search ex- 
periment performed at Kamioka Observatory with new limits in the ap — an 
(WIMPs-nucleon couphngs) plane. 



2 Experimental set-up and measurements 

The laboratory for this dark matter search experiment was newly caved near 
the Super-Kamiokande detector at Mozumi Mine of Kamioka Mining and 
Smelting Co. located at Kamioka-cho, Gifu, Japan. The laboratory is at the 
depths of 2700 m.w.e and the muon flux is reduced about five orders of mag- 
nitude compared to the surface laboratory. As the radon concentration in 
the mine air is high and has a seasonal change (40 (winter) ~ 1200 (summer) 
Bq/m^), the fresh air from the outside of the mine is blown into the labora- 
tory with a rate of 1 m^/min. The radon concentration in the laboratory is 
less than 50 Bq/m^. 
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The detector used in this dark matter search experiment is a LiF bolome- 
ter array consisting of eight LiF bolometers with a total mass of 168g. Each 
bolometer consists of a 2x2x2 cm^ LiF crystal and a neutron transmutation 
doped (NTD) Ge thermistor (1.5x1x1 mm=^ )[7] glued to the LiF with GE 
varnish. We refer to each bolometer as Dl, D2,..., and D8. Four of the LiF 
crystals are purchased from BICRON Co. Ltd. (expressed as BIG in Table 
1) while the others are from OHYO KOKEN KOGYO Co. Ltd. (expressed 
as OKEN in Table 1). NTD3 thermistors, which have been used for most 
of the previous measurements, are used for six out of the eight bolometers, 
while NTD2 and NTD4 are used for one bolometer each. Characteristics of 
the NTD thermistors are shown in Ref. [7]. Radioactive contamination in the 
crystals were previously checked with a low-background Ge spectrometer and 
was found to be less than 0.2 ppb for U, 2 ppb for Th, and 20 ppm for K. Sur- 
faces of the LiF crystals were etched with perchloric acid in order to eliminate 
remaining grains of the polishing powder. Features of the eight bolometers are 
summarized in Table 1. 

Schematic drawing of the bolometer array is shown in Fig. 1. The array consists 
of two stages, each stage holding four bolometers. The bolometers are fixed 
with spring pins and Delrin balls. Thermal contacts between the LiF crystals 
and the oxygen free high conductivity (OFHC) copper holder are achieved by 
annealed OFHC copper ribbons of 10 mm width and 0.1 mm thickness. The 
etching process and the assembly were performed in a clean bench of class 100. 
The bolometer array is encapsulated in a low-radioactivity lead (older than 
200 years) of 2cm thickness. The OFHC copper holder and the inner shield 
were etched with nitric acid in order to eliminate the surface contamination. 



Detector 


Dl 


D2 


D3 


D4 


D5 


D6 


D7 


D8 


mass (g) 


20.98 


20.98 


20.99 


20.92 


21.09 


21.07 


21.01 


21.00 


LiF manufacturer 


BIC 


BIC 


BIC 


BIC 


OKEN 


OKEN 


OKEN 


OKEN 


thermistor 


NTD2 


NTD3 


NTD3 


NTD3 


NTD4 


NTD3 


NTD3 


NTD3 


Etching (/im) 


21 


22 


26 


29 


22 


24 


27 


30 



Table 1 



LiF bolometers used for this underground measurement. 



The bolometer array and the inner shields are suspended with three Kevlar 
cords from the bottom of the dilution refrigerator which is cooled down to be- 
low lOmK. Thermal connection to the dilution refrigerator is achieved by four 
annealed OFHC copper ribbons of 10 mm width, 120 mm length, and 0.1 mm 
thickness. Each Kevlar cord is 8 cm long and eliminates the microphonic noise 
due to the helium liquefier. The dilution refrigerator is constructed with the 
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low-radioactivity materials assayed with a low-background Ge spectrometer 
in advance. 

The refrigerator is shielded with 10 cm of OFHC copper, 15cm of lead, Ig/cm^ 
of boric acid sheet, and 20 cm of polyethylene as shown in Fig. 1. Materials for 
the passive shields were previously radio-assayed with a low-background Ge 
spectrometer [8]. An air tight plastic bag is set between the copper shield and 
the cryostat. The plastic bag is filled with nitrogen gas evaporated from liquid 
nitrogen. The radon concentration was less than 0.2 Bq/m^ throughout the 
measurement period, which makes negligible background to this dark matter 
search experiment. 

Each thermistor is DC biased through a 300 MQ load resistor. The bias cur- 
rent of each bolometer is listed in Table 2. The voltage change across the 
thermistor is fed to the voltage sensitive amplifier placed at the 4K cold stage. 
A junction field effect transistor (J-FET), NEC 2SK163, is used for the cold 
stage amplifier with a gain of 10. The signal from the cold stage amplifier is in 
turn amplified by the voltage amplifier placed at the room temperature with 
a gain of 4.86 xlO^. The signal from the voltage amplifier undergoes through 
a low pass filter with a cut-off frequency of 11 Hz and is digitized with a sam- 
pling rate of 1 kHz. All the data are recorded continuously without any online 
trigger for a complete offline analysis. The data size for an hour measurement 
is 74M bytes. 

The measurement was performed from November 22, 2001 through January 
12, 2002. Eight of the eight bolometers worked, however, Dl and D5 showed 
poor signal to noise ratios because of the thermistors (NTD2 and NTD4). We 
analyzed the data obtained with the other six bolometers. 

We applied an offline level trigger to the data averaged over ±10 msec. Trigger 
levels are determined so that the trigger rates do not exceed 1 Hz. The trigger 
levels and the trigger rates are shown in Table 2. We determined the energy 
of each event by fitting the scaled template signal to the raw data. The elec- 
tric noise events and the events that are thought to have occurred within the 
thermistors are eliminated by the wave form analysis. Then the events above 
the threshold energy that is determined as the energy with efficiency=0.8 are 
selected as real bolometer signals. The trigger and event selection efficiencies 
are estimated by adding the scaled template signals to the raw data as de- 
scribed in Ref. [6]. Finally, the real bolometer events in coincidence with two or 
more detectors within ±10 msec are discarded since WIMPs hardly interact 
in two or more detectors simultaneously. Energy resolutions of the bolometers 
are estimated by the same method as the efficiency estimation. The threshold 
energies used for the event selection and estimated energy resolution near the 
threshold are listed in Table 2. 
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The calibrations arc carried out by irradiating ^°Co and ^^^Cs gamma-ray 
sources suspended between the heavy shield and the cryogenics. Compton 
edges are used and good agreements with the simulations are seen. Good lin- 
earities are seen as shown in Ref . [6] . Linearities down to 60keV are confirmed 
previously [5]. Recorded pulse heights to a given signal for six detectors are 
shown as 'gain' in Table 2. Gains are within the range of factor 7 which seems 
to be reasonable, considering that the assembly conditions are not exactly the 
same for the six bolometers. Stability of the detector response and linearities 
are checked every three or four days and the gains are found to be stable 
within ±5%. 



Detector 


D2 


D3 


D4 


D6 


D7 


D8 


Bias Current [n A] 


0.05 


0.04 


0.05 


0.05 


0.09 


0.07 


Gain[mV/keV] 


0.62 


1.1 


2.1 


1.2 


0.38 


2.5 


Trigger level [mV] 


8 


8 


15 


6 


3 


5.5 


Trigger rate [Hz] 


0.15 


0.12 


0.57 


0.011 


0.35 


0.14 


Threshold [keV] 


44.1 


24.1 


60.8 


15.7 


46.8 


8.0 


Event rate [10^^ Hz] 


0.71 


1.2 


1.0 


0.77 


0.51 


1.4 


Resolution (FWHM)[keV] 


13 


8.0 


10 


3.5 


7.8 


1.9 


Live Time [days] 


33.2 


33.4 


30.1 


34.3 


31.7 


33.3 



Table 2 

Detector responses of the six bolometers. 



3 Results 

The spectra obtained with the six bolometers are shown in Fig. 2. One of the 
spectra (D3) obtained in the pilot run at Nokogiriyama underground cell[6] is 
shown for comparison. Live times of the six bolometers are listed in Table 2. 
The total exposure is 4.1 kg-days. Background levels of D6, D7, and D8 have 
decreased by nearly one order of magnitude compared to the pilot run. D4 
suffers from electronics noises, therefore the background level of D4 is higher 
than those of other detectors. 

Prom the obtained spectra, we derived the spin-dependent WIMP-proton elas- 
tic scattering cross section (o'wiMP-p) limits. The limits are shown in Pig. 3. We 
calculated the limits in the same manner as used in Ref. [6] and Ref. [10]. Here 
we used the astrophysical and nuclear parameters hsted in Table 3. We used 
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the earth velocity which corresponds to the velocity in the measurement pe- 
riod. We used the \^J{J + 1) values calculated by assuming the odd group 
shell model[ll]. For the WIMPs heavier than 50 GeV, the limits have been 
improved by factor 5 from the pilot run. Since we determined the energy 
threshold (efficiency=0.8) higher than the pilot run, we did not improve the 
cross section limits for the lighter WIMPs compared to the pilot run. The best 
hmit is 23pb for Mwimp =40GeV. 



Dark matter density p-Q = 0.3GeVc~^cm~^ 

Velocity distribution Maxwellian 

Maxwellian velocity parameter i;o=220kms~^ 

Escape velocity v^sc =650kms~-'^ 

Earth velocity i;E=217kms~^ 

spin factor(^Li) AV(J + 1)=0.411 

spin factor(i9F) AV(J + 1)=0.647 

Table 3 

Astrophysical and nuclear parameters used to derive the exclusion limit. 



We also derived the spin-dependent WIMP-proton and WIMP-neutron elastic 
scattering cross section limits from ^Li and ^^F independently in the same 
manner as described in Ref. [9], or using 

^SDlim _ ^SDlim PWIMP-p / '-'p(N) 

<'^WIMP-p(N) — <^WIMP-N , 2 /^SD 

PWIMP-N \ '-'p 




and 



2 /r'SD \ -1 

^SDlim _ ^SDlim A^WIMP-n / ^n(N) \ /„x 

<^WIMP-n(N) — <^WIMP-N~2 7^SD~ j ' ^ > 

A^WIMP-N \ *-^n / 



where o"^^,5^_p(N-| and o"wiMP-n{N) ^^e spin-dependent WIMP-proton and 
WIMP-neutron cross section limits obtained on the assumption that all events 
are due to WIMP-proton and WIMP-neutron elastic scatterings in the nucleus 
N, respectively, cr|(^^_]vj is the WIMP-nucleus cross section limit which is di- 
rectly set by the experiments, /xwimp-p and /^wiMP-n are the WIMP-proton 
and WIMP-neutron reduced masses, respectively, /U-wimp-n is the WIMP- 
nucleus reduced mass, C'p|^^ and C^^^ arc the proton and neutron contri- 
butions to the total enhancement factor of nucleus N, respectively, and 
and are the enhancement factors of proton and neutron themselves, re- 
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spectively. Ratios of the enhancement factors calculated by assuming the odd 
group shell model are shown in Table 4. 

The enhancement factor and the cross section are related by the following 
equation: 



"^WIMP-N — ^Cf/ZwIMP-nC'n^ 

= 4GfA''Wimp-n-:(«p('S'p(n)) + an(<S'n(N)))^ — J—-, (3) 

where cr|^]y[p_j^ is the WIMP-nucleus cross section, is the Fermi couphng 
constant, is total enhancement factor of nucleus N, and the 
WIMP-proton and WIMP-neutron couplings, respectively, (5'p(N)) and (5'n(N)) 
are the expectation values of the proton and neutron spins within the nucleus 
N, respectively, and J is the total nuclear spin. (5'p(N)), (>S'n(N)), and J can be 
calculated with shell models and are shown in Table 4. 

From the obtained WIMP-proton and WIMP-neutron elastic scattering cross 
section limits, we derived the limits in the Op — an plane. We again followed 
Ref. [9]. The allowed region in the Cp — plane is defined as follows: 

\ ^ 

Op On \ TT 



/_SDlim /^SDlim I 24G^U?,.tajo 

.V^WIMP-pCN) V WIMP-n(N)/ ^^""^F/^WIMP 



< 777772-2 ' (4) 



where the relative sign inside the square is determined by the sign of (-S'p(N))/ ('S'n(N))- 

Obtained limits in the Op — plane for WIMPs with mass of 50 GeV are 
shown in Fig.4 (a). The regions between the two parallel lines are the allowed 
regions from ^Li and ^^F alone, respectively. It is seen that a single nuclide 
is not sufficient to set a finite limit in the Cp — On plane. Since we have ^Li 
and ^^F as the targets, we can take the intersection of the two regions as 
the combined allowed region from the LiF detector. This is the merit of the 



Isotopes 


J 


('S'p(N)) 


('S'n(N)) 


S(N)/*^P 


'-"n(N)/'-^n 




3/2 


0.497 


0.004 


5.49 X 10-^ 


3.56 X 10-^ 


19p 


1/2 


0.441 


-0.109 


7.78 X 10-1 


4.75 X 10-2 


23Na 


3/2 


0.248 


0.020 


1.37 X 10-1 


8.89 X 10-4 


127j 


5/2 


0.309 


0.075 


1.78 X 10-1 


1.05 X 10-2 



Table 4 

Values of (S'p(N)), (5n(N)), C'p(N)/C'p°, C'n(N)/^n° varlous isotopes. Values for '^lA 
and i^F are taken from Ref. [17] and those for ^^Na and i^^j are taken from Ref. [18] 
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detector consisting of more than two nuclides with sufficient spins. Wc show 
the aUowed regions in the ap — plane obtained by this experiment for light 
WIMPs (Fig.4 (b)) and heavy WIMPs (Fig.4 (c)). The hmits set by this 
experiment for MwiMP=30GeV are |ap| < 32 and |ap| < 133. 

Wc compared our results with those of other experiments. As the DAMA group 
had not shown the results from ^^Na and ™\ indcpcndcntly[l,13,4], we were 
not able to compare our results with the spin-dependent annual modulation 
allowed regions which is shown in Fig. 3 [4]. Instead, we compared our results 
with those of the UKDMC experiments[10,14] . Since UKDMC shows the hmits 
from ^^Na and ^^^I independently [10], we can calculate the limits in the Op — 
an plane from ^^Na and ^^"^I and combine them without any assumptions. 
It should be noted that limits in the Op — plane from the Nal detector 
are largely affected by the sensitivity ratio of ^^Na and ^^^I; the threshold 
of the detector and the quenching factor of ^^^I largely change the hmits. 
Furthermore, accurate calculation of the annual modulation allowed region 
requires all of the obtained experimental data with time information, which 
we do not have. Therefore, we did not take the risk of deriving the limits or 
allowed region in the ap — an plane from the experiments of which independent 
results from ^^Na and ^^^I are not shown. 

In Fig. 5, the allowed regions for the (a): sufficiently light WIMPs for which ^^'^I 
loses the sensitivity and Nal effectively becomes a detector with one nuclide, 
(b): WIMPs with mass near the center value of DAMA's annual modulation 
signal, and (c): WIMPs heavy enough for which ^^Na and ^^^I both have 
good sensitivities, are shown. In the three plots, the limits from Na and I 
of UKDMC[10] are shown in dotted and dash-dotted lines, respectively, Nal 
combined limits of UKDMC [10] arc shown in thin solid lines and LiF combined 
limits from this work are shown in thick sohd lines. UKDMC published the 
recalculated data[14] in which only Nal combined limits are shown, and we 
derived the limits in the ap — an plane on the assumption that the sensitivity 
ratios of ^^Na and ^^^I are the same as those shown in Ref. [10]. These new 
limits from Ref. [14] are shown in dashed lines just for reference. 

It is seen that our results can set limits complementary to the results obtained 
with Nal detectors. For the light WIMPs (lighter than 20 GeV for most of the 
detectors) to which ^^''I loses sensitivity, this is the only experiment that can 
set a finite allowed region in the ap — an plane. 

For the WIMPs with a mass of 50 GeV, which is near the center value of 
DAMA's annual modulation signal, we exclude more than two thirds of the 
parameter space allowed by the UKDMC experiments. As the limits from 
the UKDMC experiments are close to the DAMA's spin-dependent allowed 
region, as shown in Fig. 3, we should also be able to compare our limits in the 
ap — an plane. We, however, have to wait until DAMA pubhshes relevant data 
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to calculate the allowed region in the Op — On plane. 

For the WIMPs with a mass of 100 GeV, we only exclude a part of the UKDMC 
allowed region. However, it must be stressed that though our crYj^]yjp_p limits 
are about two orders of magnitude worse than the hmits of UKDMC as shown 
in Fig. 3 , we are able to tighten their limits in the a-^—a^ plane. This is because 
the spin factors of ^^F are large and the (5'p(N))/('S'n(N)) sign of ^^F is different 
from those of ''Li, ^^Na, and ^^''I. 



4 Discussions and prospects 

Though we are able to tighten the limits from the UKDMC experiment in the 
Op — an plane, the background is still high for the direct detection of WIMPs. 
We studied the background sources and found that potassium contamination 
on the surface of the LiF crystals, tritium produced by the neutron capture 
of ^Li during the development in the surface laboratory, and the uranium and 
thorium contamination in the OFHC copper holder and in the inner lead shield 
dominate the background in the energy region between 10 keV and 90 keV. 
We are planning to etch the LiF crystals more deeply in order to eliminate 
the contamination on the surface of the LiF crystals. We also plan to use NaF 
crystals as a bolometer in order to eliminate the effects of tritium. It is also 
expected that NaF crystals give more stringent limits in the Op — On plane 
since ^^F and ^^Na set limits orthogonal to each other in the Op — On plane as 
shown in Fig. 4 and Fig. 5. We are planning to install active shields between 
the LiF crystals and the OFHC copper holders to reduce the background from 
the OFHC copper holder and the inner shield. With these improvements, we 
expect to improve the sensitivity of our detector by more than one order of 
magnitude. 



5 Conclusions 

We have performed the first underground dark matter experiment at Kamioka 
Observatory. With a total exposure of 4.1 kg-days, we have obtained o"|^j^p_p 
limits factor 5 improved from the pilot run at Nokogiriyama underground cell. 
We calculated the limits in the Op — plane and showed that our results 
tightened the limits from the UKDMC experiment. For the WIMPs with a 
mass of 50 GeV, we have excluded more than two thirds of the parameter space 
in the — an plane allowed by the UKDMC experiment. As DAMA's spin- 
dependent annual modulation allowed region in Fig. 3 is close to UKDMC's 
limits, we should also be able to compare our limits in the Op — plane 
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with DAMA's spin-dcpcndcnt annual modulation allowed region when DAMA 
publishes relevant data to calculate the allowed region in the Op — On plane. 
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Fig. 1. Schematic drawings of the bolometer array. 
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Fig. 2. Low energy spectra obtained with the six bolometers are shown in thick 
soUd lines. One of the spectra (D3) obtained in the pilot run at Nokogiriyama 
underground cell [6] is shown in dashed lines. 
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Fig. 3. LiF combined 90% C.L. CyJ^j^p^p limits as a function of Mwimp are shown 
in a thick solid line. Limits derived from Li and F alone are shown in dashed lines. 
(Limits from F are actually not seen because they almost coincide with the LiF 
combined limits shown in the thick solid line.) Limits from the pilot run and from 
other experiments[13,14,15,16] are shown in a thin solid and dotted lines, respec- 
tively. DAMA's allowed region investigated in a mixed coupling framework is shown 
as "case c" in Ref. [4] is shown in a thin solid line. The MSSM predictions are also 
shown. 
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Fig. 4. Limits in the Op — On plane.: (a) For Mwimp = 50GeV, limits set by Li(F) 
alone are shown in dashed(dash-dotted) lines while the combined limits are shown 
in solid lines, (b) Limits for the light WIMPs. (c) Limits for the heavy WIMPs. 




Fig. 5. Comparison of the limits in the Op — an plane with the UKDMC data[10] for 
various WIMPs masses. 



